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ABSTRACT 

The widespread detection of synchrotron X-ray emission from the jets of low-power, nearby 
radio galaxies implies the presence of electrons at and above TeV energies. In this paper 
we explore the possibility that the TeV 7-rays detected from the radio galaxies Cen A and 
M87, which both have bright, well-studied X-ray jets, are produced at least in part by inverse- 
Compton scattering of various photon fields by the high-energy electrons responsible for the 
synchrotron X-rays on kiloparsec scales. We describe a new numerical code that we have 
developed to carry out inverse-Compton calculations taking account of all the relevant physics 
and using detailed models of the jets and the photon fields in which they are embedded, and 
show that existing constraints on the very high-energy (VHE) 7-ray fluxes of these two objects 
already place significant constraints on the magnetic field strengths in the jet in Cen A. Finally, 
we discuss the prospects for constraints on radio galaxy jet physics that may be obtained from 
observations with the Cerenkov Telescope Array (CTA). 
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1 INTRODUCTION 

The existence of X-ray counterparts to the bright radio jets in the 
nearby radio galaxies Centaurus A and M87 has been known for 
many years (Feigelson et al. 1981; Biretta, Stern & Harris 1991). 
More recently, it has become clear that such X-ray jets are common, 
possibly ubiquitous (Hardcastle et al.2001; Worrall et al. 2001) in 
classical low-power, Fanaroff & Riley (1974) class I (FRI) radio 
galaxies. From the connection to the radio through infrared, op- 
tical and ultraviolet spectrum, it appears that the X-ray emission 
is synchrotron in origin (Hardcastle et al.2001; Hardcastle, Kraft 
& Worrall 2006); for plausible (equipartition) values of the mag- 
netic field strength in the jet, this means that we are seeing emis- 
sion from relativistic electrons with TeV energies. Because of the 
synchrotron loss timescales of these electrons in the equipartition 
magnetic fields, which are of the order of tens to hundreds of years, 
it is widely accepted that the X-ray detections of well-resolved, 
kpc-scale emission from FRI jets requires an in situ acceleration 
process, and observations of Cen A and M87, the best-resolved and 
best-studied of the X-ray jets, have been used to argue both for lo- 
calized particle acceleration at shocks and for a more distributed ac- 
celeration process (Hardcastle et al. 2003; Perlman & Wilson 2005; 
Goodgeretal. 2010). 

However, these arguments are based on assumptions about the 
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magnetic field strengths in the jets. Another consequence of the 
availability of sensitive X-ray observations of radio galaxies has 
been the widespread detection of X-ray inverse-Compton emission 
from the lobes and hotspots of the more powerful FRII radio galax- 
ies (e.g. Harris et al. 1994; Hardcastle et al. 2002; Isobe et al. 2002; 
Hardcastle et al. 2004; Kataoka & Stawarz 2005; Croston et al. 
2005) which has allowed us to show that the field strengths in those 
regions of radio galaxies are generally within a factor of a few of the 
conventional equipartition value, assuming no energetically domi- 
nant proton population. But in FRI jets, the strong synchrotron X- 
ray emission means that we cannot use inverse-Compton X-rays 
to measure magnetic field strength; indeed, with the exception of 
the recent Fermi measurement of magnetic field strengths in the 
giant lobes of Cen A (Abdo et al. 2010) we have no direct inverse- 
Compton constraints on the magnetic field strength of any compo- 
nents of the FRI population. Among other things, this restricts our 
ability to draw conclusions about the particle acceleration proper- 
ties of FRI jets. 

The inference from the X-ray observations that very-high- 
energy leptons are present in the FRI jets does, however, imply that 
inverse-Compton emission from these regions should be detectable 
up to very high energies: electrons with energies in the 1-10 TeV 
energy range can inverse-Compton scatter a suitable photon popu- 
lation into the TeV 7-ray band. FRI jets are in principle well sup- 
plied with parent photon populations. In addition to the cosmic mi- 
crowave background, which is the principal photon field involved 
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in the inverse-Compton detection of radio galaxy lobes, and the ex- 
tragalactic background light (EBL), which is also always present, 
synchrotron self-Compton emission (SSC) may be important due 
to the relative compactness and brightness of FRI jets, and they are 
also exposed to any emission that may be visible to them (though 
usually not directly to us) from the active nucleus itself and from 
the parsec-scale jet, which will normally be beamed in the direction 
of the kpc-scale jet. Finally, in the centres of elliptical galaxies, the 
energetically dominant photons are starlight: the energy density in 
starlight at the centre of a large elliptical galaxy can be ~ 2 x 1CP 12 
J m -3 , exceeding that in the z — CMB by nearly two orders of 
magnitude. This fact led Stawarz, Sikora & Ostrowski (2003) to 
predict that the inverse-Compton scattering of starlight would pro- 
vide an important contribution to the emission of FRI jets in the 
TeV regime. 

The opportunity to confront these inverse-Compton predic- 
tions has been provided by the detection of M87 and, recently, Cen 
A with Cerenkov imaging telescopes sensitive in the TeV range 
(Aharonian et al. 2003; Aharonian et al. 2009). The spatial resolu- 
tion of the current generation of these instruments is not sufficient 
to distinguish between emission from the kpc-scale jet and emis- 
sion from the active nucleus or from the jet on smaller scales; from 
the widespread detection of, often strongly variable, TeV emis- 
sion from blazars we know that small-scale jets are certainly capa- 
ble of producing high-energy 7-rays, while the detection of short- 
timescale variability in the TeV emission from M87 (e.g. Acciari et 
al. 2009, 2010) is evidence that some at least of this emission orig- 
inates on small spatial scales, either in the pc-scale jet or possibly 
in the peculiar jet knot HST-1 (Cheung, Harris & Stawarz 2007). 
However, even if we treat the level of TeV emission (or, in the case 
of M87, its non- variable component) as giving us an upper limit on 
the inverse-Compton emission from the kpc-scale jets, these detec- 
tions give us a lower limit on the magnetic field strength that cannot 
be obtained in any other way, so long as we can make an accurate 
model of the inverse-Compton emissivity in the TeV regime. 

Detailed inverse-Compton modelling of kpc-scale jets is chal- 
lenging because of the large number of photon fields that need to be 
taken into account (as discussed above) and because of the strong 
point-to-point variation in the photon energy density in several of 
them; it is clear that the traditional methods of constructing one- 
zone synchrotron/inverse-Compton models are unlikely to give an- 
swers that can be relied on to better than an order of magnitude. 
In this paper we present a framework for inverse-Compton mod- 
elling of M87 and Cen A and show that the existing TeV detec- 
tions already provide quite strong constraints on the magnetic field 
strength in these objects. We also discuss the prospects for better 
constraints on these objects, and perhaps detections of others, with 
the next-generation capabilities to be provided by the Cerenkov 
Telescope Array (CTA). 

The remainder of the paper is organized as follows. In Section 
[2] we provide an overview of the code that we use to carry out the 
jet modelling and inverse-Compton calculations, and in Section [3] 
we comment on some of the limitations imposed by our approach. 
Section [4] describes the general approach we take to matching our 
jet models to the data and the specific constraints that we use for 
Cen A and M87. Section[5]gives the results of our modelling, Sec- 
tion [6] describes our view of the outlook for future facilities, and 
conclusions are given in Section|7] 

Throughout the paper we assume a distance of 3.7 Mpc to Cen 
A (the average of 5 distance indicators presented by Ferrarese et 
al. 2007; a slightly higher distance, 3.8 Mpc, is suggested by later 
work, e.g. Harris et al. 2010, but we retain the lower distance for 



consistency with our earlier work, noting that it makes no signifi- 
cant difference to our results) and a distance of 16.4 Mpc to M87 (a 
weighted mean of 4 estimators presented by Bird et al. 2010). For 
objects at larger distances we assume H = 70 km s _1 Mpc -1 , 
Q M = 0.3, and Ov = 0.7. 



2 THE CODE 

The code we use is based on the framework described by Hardcas- 
tle et al. (2002: hereafter H02). In that paper we were concerned 
primarily with the X-ray synchrotron-self-Compton emission from 
the bright, spatially resolved hotspots of FRII radio galaxies: the 
new feature of the code described there was the fact that it allowed 
us to consider any spatial distribution of the scattering electrons, 
rather than, as previously, restricting ourselves to uniformly filled 
spherical regions. Such an approach is also required for our current 
problem since, as discussed above, the photon fields illuminating 
the jet are strongly dependent on position, so that one-zone mod- 
els will be completely inadequate. Because the use of this code is 
the main new feature of the paper, in this section we describe it in 
detail. 

The basic principle of all parts of the code is the idea that 
spatially resolved structures can be modelled in terms of three- 
dimensional Cartesian 'grids' of finite-sized, cuboidal volume el- 
ements, with each volume element constituting a separate 'zone' 
of the synchrotron and inverse-Compton modelling and having one 
or more rest-frame physical values associated with it. So, for ex- 
ample, the synchrotron emission from the jet is modelled by two 
such grids, one containing values for the electron energy spectrum 
normalization and one for the magnetic field strength. (A third grid, 
used by H02 but not in the present work for reasons that will be dis- 
cussed below, specifies which of a finite number of electron spectral 
models describe each region.) Computing a quantity such as the to- 
tal synchrotron flux density from the source then involves doing 
a separate synchrotron emission calculation for each non-zero el- 
ement of the grid and generating another grid, in this case one of 
synchrotron volume emissivity at a specified frequency, which can 
then be summed over to give the luminosity or flux density from the 
modelled object (transformed into the observer's frame). In addi- 
tion, importantly, the synchrotron emissivity grid can be projected 
on the plane of the sky to give an image that is proportional to the 
observed synchrotron surface brightness. The main operations of 
the code can all be described in terms of operating on a parameter 
file and, ordinarily, one or more grids and generating a new grid, 
which may then be operated on in turn. 

This structure means that, in order to compute inverse- 
Compton properties of the modelled source, we need to compute 
the inverse-Compton emissivity for each element of the grid sep- 
arately. Doing this, of course, is what gives us the ability to make 
models of the source with adequate resolution to represent its real 
physical behaviour, but it comes at the cost of greatly increased 
computational requirements, since to achieve adequate resolution 
(which in practice means having a cell size significantly smaller 
than the jet radius and much smaller than the characteristic scale 
of the stellar emissivity profile) our grid needs to have ~ 10 4 non- 
zero volume elements distributed throughout the jet. The inverse- 
Compton kernel that we use is the one presented in great detail 
by Brunetti (2000): while in H02 we used the approximations 
given by Brunetti appropriate for Thomson scattering in the rest 
frame, in this paper we use the full version with no approxima- 
tions, so as to take account of Klein-Nishina effects, which adds 
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to the computational complexity of the problem. Specifically, our 
inverse-Compton emissivity, assuming an isotropic electron pop- 
ulation within each cell, is given by Brunetti's equations 31 and 
32, and we find the minimum Lorentz factor capable of scatter- 
ing between two given energies using equations 33 and 13. The 
emissivity is then integrated over all available photon and elec- 
tron energies. Since jets are beamed, we also need to take account 
of the effects of special relativity in the inverse-Compton calcula- 
tions. For lab-frame isotropic photon fields such as the CMB or 
the EBL, we can simply integrate over the whole sky, taking ac- 
count of the Doppler shift and transforming all angles into the 
jet rest frame, to obtain a conversion between the normalization 
of the electron energy spectrum in a given cell and the inverse- 
Compton volume emissivity. However, in the cases of SSC and 
inverse-Compton scattering of starlight, where the photon field is 
anisotropic, the problem is harder: we need to compute the illumi- 
nation of every grid element by a large number of other grid ele- 
ments, and in principle each one has a different Doppler factor and 
a different jet-frame inverse-Compton scattering angle. We use a 
modified version of the simplification described by H02, in which 
we first calculate the integral of Brunetti's eq. 31 over a reasonably 
well-sampled lookup table of Doppler factors and scattering angles 
and then interpolate over this lookup table to find the contribution 
to the inverse-Compton emissivity for every illuminating cell and 
for every illuminated cell. While this process is still computation- 
ally expensive, it is very much less so than it would be if we were 
to do the full integral for every cell, although it does require the 
use of some simplifying assumptions, which are described in more 
detail below. The SSC and starlight-illumination codes make use of 
the MPI framework, as implemented in MPICH23. to allow the cal- 
culation of the lookup table and the actual emissivity computation 
to be distributed over a large number of machines: while individual 
calculations are practical on modern multi-core desktop machines, 
the final computations in which jet speed and angle to the line of 
sight were allowed to vary were done on the University of Hert- 
fordshire clustefl typically using 16 physical compute servers and 
128 cores. 



3 LIMITATIONS AND APPROXIMATIONS 

The models we use have the following restrictions and limitations. 

(i) We assume a single electron energy spectrum throughout 
the jets. Allowing for a discrete number of different spectra, as in 
H02, or, worse, a continuously spatially varying electron spectrum 
throughout the jet would prevent us from applying the computa- 
tional simplification described above - in effect, we would have to 
do the full numerical integration for the illumination of every cell 
by every other cell in the SSC and starlight cases. At present this 
is computationally intractable. Since we know that the synchrotron 
spectra of the jets actually do vary as a function of position, at least 
at high energies (e.g. Hardcastle et al. 2007) this is not an ideal 
approximation, but it is necessary to make progress. 

(ii) We assume a single value of the magnetic field strength 
throughout the jet. This is also needed to simplify the synchrotron 
self-Compton calculations; without it, the illuminating synchrotron 
spectrum from each region of the jet would be different, so again 
we would greatly increase the required number of integrals over the 

1 See http://www.mcs. anl.gov/research/projects/mpich2/ 

2 http://star. herts . ac . uk/progs/computing .html 



Brunetti equations. This constraint coupled with the previous one 
means (a) that the only means we have of modelling the surface 
brightness variation of the jet is to vary the normalization of the 
electron energy spectrum, and (b) means that we make the implicit 
assumption that the jet looks the same at all frequencies - which 
again we know to be an approximation that breaks down at the 
highest energies. Because the equipartition field strength depends 
only weakly on electron energy density, this approximation is not 
likely to be seriously problematic for realistic jet models. 

(iii) We assume a uniform velocity field throughout the jet - all 
components of the jet are moving at the same speed and the same 
direction. We know from modelling such as that of Laing & Bridle 
(2002) that this is not likely to be accurate, and that the magnitude 
and direction of the velocity field are likely to vary radially and 
along the jet. In fact, incorporating a varying velocity field into 
our modelling would not be too difficult, but we have little or no 
constraint on anything but the bulk jet speeds for our targets, and 
so have not implemented this feature. The effects of beaming on our 
results are in fact relatively modest, given the jets' low speeds and 
comparatively large angles to the line of sight (see discussion of 
M87 below, Section [5^2l . and so we do not expect this assumption 
to have any very important effects. 

(iv) We assume that the starlight volume emissivity distribution 
is spherically symmetrical and can be modelled by a Mellier & 
Mathez (1987) distribution, based on the deprojection of a de Vau- 
couleurs (r 1//4 ) profile. More importantly, we assume, for the same 
reasons as for the previous point, that the emission spectrum from 
each region of the host galaxy is the same. In the case of M87, 
this is not too bad an approximation; in the case of Cen A, we 
know that the dust emissivity must in reality be very differently 
distributed from that of the starlight. However, both are strongly 
centrally peaked, which is the most important factor from the point 
of view of their effect and the jet, and we have not considered it 
necessary, for example, to make a two-component model of the dust 
and starlight emissivity in Cen A, although the structure of the code 
means that that would be possible in future if a more detailed model 
were warranted. 

(v) We neglect any illumination from the hidden AGN or blazar. 
Self-illumination of the modelled part of the jet is obviously taken 
account of by the SSC modelling, but we have very little informa- 
tion on the structure and intrinsic luminosity of the pc-scale jet and 
even less on the spectrum of the AGN as seen by the kpc-scale 
jet. This in principle means that any flux/luminosity calculations 
we consider should probably be taken to be lower limits. In prac- 
tice we may not be making too large an error by neglecting these 
photon fields: both are Doppler-suppressed in the frame of the kpc- 
scale jet and both suffer strongly from inverse-square dilution with 
distance from the centre of the galaxy, which is not the case for the 
starlight or SSC fields. 

(vi) We neglect photons from the diffuse, kpc-scale lobes of our 
targets: these are not likely to be energetically important compared 
to the other photon fields and the comparatively low-frequency 
photons they predominantly emit will be hard to scatter to TeV en- 
ergies. 

(vii) Finally, we do not attempt to model the attenuation of TeV 
7-rays by photon-photon interactions within the host galaxy (our 
objects of interest are close enough that interactions with the EBL 
as the 7-rays travel between the sources and the detector may safely 
be neglected). The cross-section for interactions of TeV photons 
peaks at target photon energies corresponding to the near-IR, where 
the photon density is relatively low, and modelling has shown that 
attenuation by host-galaxy starlight can be safely neglected both for 
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the Milky Way (Moskalenko, Porter & Strong 2006) and for more 
strongly star-forming galaxies (Gilmore & Ramirez-Ruiz 2010). In 
the specific case of Cen A, Stawarz et al. (2006a) have shown that 
the attenuation is at the ~ 1 per cent level. We therefore feel justi- 
fied in neglecting this effect. 



4 MODELLING THE JETS AND HOST GALAXIES 

Given the constraints above, our modelling of the jets involves 
attempting to reproduce (1) the observed overall jet synchrotron 
spectral energy distributions and (2) the approximate appearance of 
the jets in radio synchrotron emission. Because we know that high- 
energy electrons are necessary to produce high-energy photons, we 
only model in this way the region of the jets that are detected in 
X-rays. 

We know that the integrated spectra of these jets tend to 
be adequately modelled as broken power-laws (Hardcastle et al. 
2001, 2006) and so objective (1) is achieved by determining the 
low-energy electron energy index from observations in the radio 
through optical band, by taking the X-ray spectral index to indicate 
the high-energy electron energy index and by adjusting the overall 
normalization and the energy of the break, 7brcak, so as to ensure 
that the synchrotron spectrum for the jet always passes through 
the radio and X-ray data points. At the same time, we adjust the 
magnetic field strength in the jet so as to obtain overall equiparti- 
tion - that is, the total energy in electrons and magnetic field inte- 
grated over the jet are equal. This does not correspond to equipar- 
tition at every point in the jet because, as discussed above, the 
electron spectrum normalization varies from point to point in the 
grid, while there is only a single magnetic field strength through- 
out the jet. However, no point in the jet is very far from this more 
traditional equipartition condition. We assume a low-energy cutoff 
7min = 100 and a high-energy cutoff 7 max = 5 x 10 9 (i.e. corre- 
sponding to electrons with PeV energies, well above the maximum 
energy to which Cerenkov telescopes are sensitive; the results are 
not sensitive to this value as long as it is above ~ 10 TeV). 

Objective (2) is achieved by adjusting the spatial dependence 
of the electron normalization so as to give a reasonable match to 
the observed synchrotron surface brightness. We model regions of 
the jets as simple geometrical structures, such as cylinders or trun- 
cated cones, with uniform or smoothly varying electron densities 
as a function of position within the structure. Because we have to 
model the three-dimensional electron distribution, we do not at- 
tempt to reproduce the detailed structures seen in either the radio or 
the X-ray synchrotron emission, and of course (as discussed above) 
we cannot attempt to reproduce the differences between them. Our 
models are thus really intended to represent the jets rather than to 
give a detailed picture. We emphasise though that the structures that 
we do not model (principally compact knots seen predominantly in 
X-ray) are not expected to make a very strong contribution to the 
inverse-Compton flux from the sources. 

A key factor in the modelling of both jets is the jet bulk 
Lorentz factor and angle to the line of sight. Obviously these have a 
strong effect on the rest-frame energy density in both electrons and 
magnetic field, but also photons: the latter effect gives rise to the 
well-known boosting of inverse-Compton scattering of isotropic 
photon fields like the CMB in the case of highly beamed jets (e.g. 
Tavecchio et al. 2000) but the former must also be accounted for, 
in the sense that the observed radio and X-ray data points must be 
corrected for Doppler effects before the electron energy normaliza- 
tion and magnetic field strengths are corrected. Slightly less obvi- 



ously, the unknown angle to the line of sight has a strong effect on 
the volume of the jet (also affecting the magnetic field strength and 
electron energy density) and the spatial positions it occupies (which 
is important in the cases of SSC and inverse-Compton scattering of 
starlight, where the photon field depends on position). 

The major direct constraint on the bulk Lorentz factor Y and 
the angle to the line of sight 6 comes from observations of proper 
motions in the kpc-scale jets (Biretta, Zhou & Owen 1995; Hard- 
castle et al. 2003; Goodger et al. 2010). If we assume that the ob- 
served speeds correspond to bulk motions, then the apparent speeds 
/3 ap p are related to the true underlying /3 by the standard formula 

/3sin6> 

Papp — a „ 

1 — p cos 6 

So for any choice of angle to the line of sight, we can compute 
a corresponding /3 from the measured value of /3 ap p. Our initial 
approach in jet modelling was therefore to use a fixed value of /3 apP , 
determined from observations, and model the jet at all angles to 
the line of sight, computing the inverse-Compton flux density at 
1 TeV for all photon fields. This gives us a sense of the effect of 
this angle uncertainty on our TeV predictions. We then focus on 
the most likely angle to the line of sight and on more plausible 
beaming models in the later parts of the paper. 

The following subsections give details of the models applied 
to each of the two jets. 

4.1 Cen A 

The radio flux density of the Cen A jet, 30.6 Jy, is measured from 
the 1.4-GHz map of Hardcastle et al. (2006). The X-ray flux density 
at 1 keV, 134 njy, is the sum of the normalizations of the extended 
emission regions considered by Hardcastle et al. (2007). We do not 
consider the compact knots, although, as shown by Goodger et al. 
(2010), these account for about half the absorption-corrected 1-keV 
flux density of the source. In general compact features like these 
knots are rather poor sources of inverse-Compton emission where 
scattering of external photon fields is concerned, though they are 
better sources of SSC emission. In any case, we lack the detailed 
information on the knots' structure and electron energy spectrum 
needed to model them in the framework provided by our code, so 
we have chosen to model only the extended component. This means 
that our predictions for SSC emission from Cen A are probably on 
the conservative side. 

We model the electron energy spectrum using the model fitted 
to the radio, mid-IR and X-ray data by Hardcastle et al. (2006), 
which is a broken power-law in electron energy: the low-energy 
electron energy index is 2.06 steepening to 3.88 at high energies, 
corresponding to the photon index of 2.44 measured by Hardcastle 
et al. (2006) for the middle region of the X-ray jet. As described 
above, the normalization and break energy of the electron spectrum 
are adjusted to keep consistency with the observed radio and X- 
ray measurements. We now know (Hardcastle et al. 2007) that the 
photon index even of the extended component of the jet, as well 
as the radio/X-ray ratio, varies with distance along the jet, but, as 
discussed above, we need to adopt a single electron spectrum in 
order to make progress. 

We model the spatial structure of the jet as a truncated cone 
with inner radius 3 arcsec, outer radius 22 arcsec, and projected 
length 230 arcsec. Within this cone, the electron normalization de- 
creases linearly with radius and as (1 — d) where d is the distance 
along the axis of the cone. This electron distribution gives a rea- 
sonable match to the observed (knot-free) X-ray surface brightness 
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Figure 1. Comparison of the real and model synchrotron images of Cen A. Top: 5-GHz VLA map with 6-arcsec resolution showing the radio emission from 
the X-ray bright region of the jet. Middle: Chandra X-ray image from the data used by Hardcastle et al. (2007). Bottom: model synchrotron emission, assuming 
9 = 50° ; axis labels are in (projected) kpc. The model does not attempt to reproduce the knotty and patchy structure seen in the real jet (the inner part of the 
jet is dominated by knots in the radio and the X-ray) but does try to reproduce the general ran of surface brightness in the diffuse X-ray emission. The core 
(point source to the left, top and middle panels) is not represented on the model image. 



(Fig. Q3 : we constrain our models to reproduce the spatial distribu- 
tion of the diffuse X-ray synchrotron emission given that we are 
most interested in the distribution of the TeV electrons. The re- 
sults of our modelling are only weakly sensitive to our assumptions 
about the spatial distribution of the electrons. 

We take /3 apP to be 0.534, as measured by Goodger et al. 
(2010) for knot A1B, and for the host galaxy modelling we 
adopt the values used by Croston et al. (2009); flux densities for 
NGC 5128 from the ultra-violet down to the far-IR are taken from 
NED. Both for Cen A and for M87 we use the model for the EBL 
used by Hardcastle et al. (2009), i.e. the model of Raue & Mazin 
(2008). 

Finally, we take the measured VHE flux to be (1.56 ± 0.67) x 
10~ 12 photons cm~ 2 s _1 above 250 GeV, and the measured photon 
index to be 2.73±0.65, as reported by Aharonian et al. (2009). This 
means that the 1-TeV flux density of Cen A is (1.6 ± 0.7) x 10~ 16 
Jy, if we assume a fixed photon index of 2.73. 



4.2 M87 

The normalizing radio flux density of M87, 5.25 Jy, is taken from a 
5-GHz VLA radio map at 0.4-arcsec resolution. The 1-keV X-ray 
flux density we use (279 njy) and the photon index (2.21) are de- 
rived from fitting a single power-law model to 667 ks of archival 
Chandra data (Harwood & Hardcastle, in prep.). The region used 
for the spectral extraction includes the whole extended X-ray jet, 
but excludes the bright, strongly varying knot HST-1 (e.g. Harris 
et al. (2006). Clearly measuring the X-ray flux in this way is less 
conservative than our approach in the case of Cen A: in M87 we 
cannot separate discrete X-ray knots from the diffuse X-ray emis- 
sion, since the spatial resolution is so much worse. If the situation 
in Cen A were replicated in M87, then the total X-ray flux from 
the jet would exceed the X-ray flux density from the diffuse com- 
ponent only by a factor ~ 2. As with Cen A, we model the overall 
electron spectrum as a broken power law, with the same injection 
index as for Cen A and with a break that replicates the observed 
photon index in the X-ray. 

We model the jet as a combination of a truncated cone of 
length 11 arcsec, with inner radius 0.4 arcsec and outer radius 1.0 
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Figure 2. Comparison of the real and model synchrotron images of M87. Top: 5-GHz VLA map with 0.4-arcsec resolution showing the radio emission from 
the X-ray bright region of the jet. Middle: Chandra X-ray emission from the same region. Bottom: model synchrotron emission, assuming 9 = 45°; axis 
labels are in kpc. The core is not represented in the model image. 



arcsec, and a cylinder of radius 1.0 arcsec and length 6.5 arcsec, 
with electron density decreasing linearly with radius in both struc- 
tures, and with a jump in electron density of a factor 1.5 at the 
boundary between the cone and the cylinder. This roughly repro- 
duces the observed radio structure of the M87 jet (Fig. [2} with the 
boundary between conical and cylindrical structures representing 
the change in the jet opening angle and surface brightness at knot 
A. 

Jet speeds are estimated using the observations of Biretta et 
al. (1995). In this very detailed study various components of M87 
are observed to have different apparent speeds: Biretta et al. in fact 
suggest that the observations are consistent with high bulk speeds 
(7 ~ 3) in the inner jet and much slower speeds beyond knot A, 
which is consistent with the interpretation of knot A as a jet-wide 
shock. Modelling the SSC emission in this scenario would require 
us to calculate the mutual illumination of two jet regions moving 
relativistically with respect to each other, which is not yet imple- 
mented in our inverse-Compton code. For the sake of simplicity 
and for ease of comparison with Cen A we adopt /3 app = 0.479, 
which is the speed quoted by Biretta et al. for the edge of knot A. 
However, we consider below the possibility that the jet may have a 
bulk speed much higher than this. 

We model the host galaxy in the same way as that of Cen A, 
assuming that it follows a de Vaucouleurs profile, with parameters 
taken from Liu et al. (2005). For optical flux densities to determine 
the starlight SED we again take values from NED, but at IR wave- 
lengths we need to take into account the effects of contamination 
from synchrotron emission: Baes et al. (2010) have used Herschel 



data to put upper limits on the amount of cold dust that can be 
present in the host galaxy, showing that almost all the FIR emis- 
sion lies on a simple extrapolation of the known lower-frequency 
synchrotron spectrum. We base our model in the mid to far-IR on 
the modified black-body model that they fit to the residual fluxes 
after subtraction of the synchrotron component, which means that 
M87 is very weak in this region of the spectrum compared to Cen 
A. 



At TeV energies M87 is known to be strongly variable (Ac- 
ciari et al. 2009): in the past it has been argued that the flaring TeV 
emission might originate in HST-1 (e.g. Cheung et al. 2007) but 
observations of a strong flare at TeV energies that had no X-ray 
counterpart but was close to the start of a jet ejection event in the 
radio led Acciari et al. (2009) to argue that the most likely site of the 
strongly varying TeV component is in the parsec-scale jet. In any 
case, it is clear that any component of the TeV emission variable on 
timescales of days to weeks cannot be related to the diffuse emis- 
sion from the kpc-scale jet. For comparison with our model predic- 
tions we have taken the average flux density measured by Acciari 
et al. (2010) from the VERITAS observations in 2009, in which no 
significant variability of the TeV emission was observed. They re- 
port a flux of (1.59 ± 0.39) x 10" 12 photons cm -2 s" 1 above 250 
GeV with r = 2.5 (as measured for the 2008 data, which however 
include a gamma-ray flare), which is coincidentally almost identi- 
cal to what we have used for Cen A above, and which corresponds 
to a 1-TeV flux density (assuming V = 2.5) of (2.0±0.49) x 10" 16 
Jy. 
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Figure 3. Predicted 1-TeV flux density from the Cen A jet as a function of the angle made by the jet to the line of sight, 8. The three dashed horizontal lines 
show the best estimates of the 1-TeV flux density from the HESS measurements and its ltr error. 



5 RESULTS 
5.1 Cen A 

Fig.[3]shows the predicted 1-TeV inverse-Compton flux density as a 
function of jet angle to the line of sight 6 for Cen A, given the model 
discussed in Section |4~T| for an equipartition £?-field as discussed 
in Section|4] 

The general features of this plot are as expected from first prin- 
ciples. SSC emission is strongest when the jet is in the plane of the 
sky, mainly because it is physically smallest in that orientation so 
that the inverse-square suppression of SSC is least, but also because 
the correction for Doppler suppression increases both the electron 
and the photon density in the rest frame of the jet. CMB and EBL 
flux densities both increase as the angle to the line of sight gets 
smaller, both because of the increasing volume of the jet, which 
means that the total number of scattering electrons increases, and, 
at small angles, because of the effects of Doppler-boosting of these 
isotropic photon fields into the jet rest frame. The starlight com- 
ponent scales in a similar way, but increases more weakly at very 
small angles to the line of sight because for these angles the bulk 
of the jet is well away from the centre of the galaxy where the en- 
ergy density in starlight is at its peak, and because for these angles 
most of the scattering, especially when beaming effects are taken 
into account, is at small deflection angles where the cross-section 
for inverse-Compton scattering is lowest. We see that the predicted 
1-TeV flux density is very similar to what is measured using HESS 
for 20 < 6 < 30° ; at smaller angles to the line of sight we over- 
predict the TeV flux density. 



The actual angle to the line of sight for Cen A is not well 
known: Hardcastle et al. (2003) have summarized the constraints 
from radio data. If we adopt their favoured value of 50° , we see that 
the inverse-Compton prediction is 0.27 x 10~ 16 Jy, a factor ~ 6 be- 
low the flux density measured by HESS. Fig.|4]shows the predicted 
inverse-Compton SED from such a model. We see that the spec- 
tral index in VHE 7-rays is steeper than that in the X-ray, which 
is a result of the Klein-Nishina correction to the scattering cross- 
section: the predicted photon index for the total inverse-Compton 
spectrum is in fact around 3.1, but this is consistent within the er- 
rors with what is observed by HESS (note that if we assume that 
this is the correct spectral index, the 1-TeV flux density calculated 
from the HESS measurement comes down to (1.2 ± 0.5) x 10~ 1B 
Jy, a factor ~ 4 above the equipartition prediction). We conclude 
that for equipartition and for our best estimates of the jet properties 
(9 = 50°, P = 0.51) inverse-Compton scattering from the kpc- 
scale jet is unable to produce the observed TeV emission, although 
it is possible to produce all the observed 7-rays with smaller angles 
to the line of sight (we note that Hardcastle et al. 2003 suggested 
that 6 ~ 20° was required on the assumption of jet/counterjet sym- 
metry in the kpc-scale jet). 

Any departure from equipartition alters the situation signif- 
icantly. As B decreases below B cq , the electron number density 
and thus inverse-Compton emissivity increase rapidly (Fig.[5]l. The 
predicted flux goes roug hly as (B/B eq y 2A at 1 TeV: so, for the 
model with 6 = 50° discussed above, the HESS data require 
B > 0.6-Beq- In the lobes and hotspots of FRII radio galaxies (e.g. 
Hardcastle et al. 2004, Kataoka & Stawarz 2004, Croston et al. 
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Figure 4. Spectral energy distribution for the inverse-Compton emission 
for a model of the Cen A jet with 9 = 50°, /3 = 0.51 and B = B cq . The 
solid black line is the synchrotron emission, and the solid red line the total 
inverse-Compton from all modelled photon fields. The separate contribu- 
tions from the four photon fields are also plotted: from lowest to highest at 1 
TeV, these are the EBL (orange line), the CMB (blue line), SSC (green line) 
and starlight (magenta line). The red dots show our constraints on the radio 
and X-ray flux densities and the green cross is our adopted 1-TeV flux den- 
sity, with the bow-tie showing the observational uncertainties on the photon 
index. Frequencies and energies are plotted in the observer's frame. 

2005) we routinely see departures from equipartition at this level, 
so we do not regard it as at all implausible that the TeV emission 
really is wholly produced by the kpc-scale jet for 6 ~ 50°. We 
can also say that models where B is more than a factor of a few 
below £> cq are convincingly ruled out by the data, irrespective of 
8, particularly when we consider that many possible other emission 
processes may contribute to the observed TeV emission (Aharonian 
et al. 2009) and that our models of the source TeV electron popula- 
tion in Cen A are quite conservative. Stawarz, Kneiske & Kataoka 
(2006b) reached a similar conclusion based on the FRI source pop- 
ulation's contribution to the 7-ray background, but this is the first 
time that TeV emission has been used directly to constrain the mag- 
netic field strength in an individual source. 
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Figure 5. Inverse-Compton flux density from the jet of Cen A as a function 
of magnetic field strength, normalized to the equipartition value, for the 
model with 9 = 50° discussed in the text. Magnetic fields only a small 
factor below the equipartition field can allow the model to reproduce all the 
observed TeV emission. 

tition of a factor ~ 8 would be required to bring the flux up to 
observable levels. Even assuming 7 = 3 for the whole jet [as ar- 
gued by Biretta et al. (1995) for the inner jet only], while retaining 
6 = 45°, only increases the predicted 1-TeV flux density by a fac- 
tor ~ 2. We can again claim to have constrained the field strength 
from these observations, but it is a much weaker constraint, and 
our prediction is that the bulk of the TeV emission from M87 must 
come from some other source unless the angle to the line of sight 
is very small. One interesting point here is the relative weakness of 
the starlight component in M87 compared to that Cen A: although 
the energy density in starlight in M87 is significantly higher than 
that in the CMB, and this is reflected in the relative normalizations 
of the peaks in Fig. [7), the starlight contribution actually falls to 
levels comparable to that of the CMB emission by 1 TeV, whereas 
'starlight' dominates by a large factor at the same energies in Cen 
A, as can be seen in Fig. [4] We believe that the reason for this is 
the very different spectra adopted for the two host galaxies. The 
'starlight' model in Cen A includes a large contribution from dust 
emission from the central dust lane, whereas there is no evidence 
for dust emission from M87 at any significant level (Section [4. 21 ). 
Klein-Nishina effects significantly reduce the efficiency of scatter- 
ing of the predominantly optical and near-IR photons from the M87 
stellar population into the TeV band. At these energies, a simple 
comparison of the energy densities of the photon populations is not 
a reliable guide to their contribution to inverse-Compton emission: 
this also suggests that strong dust emission may be a pre-requisite 
for a bright TeV-emitting jet. 



5.2 M87 

The situation is clearly different for the M87 jet. Fig.|6]again shows 
inverse-Compton emission as a function of angle to the line of sight 
for our adopted model with /3 app = 0.479. In this model, SSC 
emission actually dominates over the other photon fields for large 
8, presumably a consequence of the fact that M87's jet is physically 
smaller than and significantly more luminous than Cen A's. The key 
difference, though, is that the net IC emission does not approach 
the levels of the 2009 VERITAS flux for 6 > 10° : for 9 = 45° (the 
SED for which is plotted in Fig. |7), it is more than two orders of 
magnitude below the observations, and a departure from equipar- 



6 PROSPECTS FOR NEXT-GENERATION FACILITIES 

The Cerenkov Telescope Array (CTA) is the next-generation TeV 
gamma-ray facility, currently in the design phase (CTA Consortium 
2010). It will both be significantly more sensitive and have sig- 
nificantly higher angular resolution than existing facilities such as 
HESS and VERITAS. Unfortunately, the likely angular resolution 
(expected to be < 1 arcmin at 1 TeV: e.g. CTA Consortium 2010) 
is too low to resolve even M87's kpc-scale jet (with a length in the 
X-ray of ~ 20 arcsec) from its active nucleus, and the same is true 
for most other kpc-scale X-ray jets at distances greater than M87's. 
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Figure 6. Predicted 1-TeV flux density from the M87 jet as a function of the angle made by the jet to the line of sight, 6. The three dashed horizontal lines 
show the best estimates of the 1-TeV flux density from the 2009 VERITAS measurements, as discussed in the text, and its la error. 



log 10 (Energy/eV) 
5 



log 10 (Energy/eV) 
5 




og 10 (Frequency/Hz) 



Figure 7. Spectral energy distribution for the inverse-Compton emission for a model of the M87 jet with B = B e q. Left panel: = 45°, (3 = 0.46. Right 
panel: 9 = 45°, 7 = 3. Lines and points as in Fig. [4] the bowtie shows V = 2.5 ± 0.3, as measured by Acciari et al. (2010) for the pre-flare data in 2008, as 
no measurement is available for the 2009 dataset. 
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Figure 8. A map of 1-TeV inverse-Compton emission from the jet of Cen A (assuming 9 = 50°, j3 = 0.51 as described in the text). The image can be seen to 
be significantly brighter close to the active nucleus and the centre of the host galaxy (left-hand side) compared to the model synchrotron image of Fig. [TJ Axis 
labels are in kpc. 



The best hope for an observational test of the idea that kpc-scale 
jets can be responsible for TeV 7-ray emission will, unsurprisingly, 
be provided by Cen A. Here the X-ray jet, which extends for ~ 4 
arcmin on the sky, should be clearly resolved by the CTA. Unfor- 
tunately in our models, in which scattering of starlight dominates 
the inverse-Compton emission, the jet is still brightest in the un- 
resolved central regions (Fig. |8}. However, if we take the model 
with 6 = 50° displayed in (Fig. [§}, and assume that all the cur- 
rently observed VHE emission is produced by the jet, then half the 
total observed 1-TeV flux density should be produced in the re- 
gion outside 1 arcmin (1 kpc), and this extended emission should 
be relatively easy for the CTA to detect, given its greatly improved 
sensitivity over its predecessors. A moderate-duration observation 
of Cen A with the CTA will therefore very probably be able to 
make a conclusive measurement of the magnetic field in an FRI jet 
for the first time. Even a non-detection would place a lower limit 
on the magnetic field strength that will be considerably higher than 
the equipartition value, and would therefore be of great interest as 
no other component of a radio galaxy for which the magnetic field 
strength has been measured exhibits this kind of behavioufl 

Other than Cen A, the best FRI jets for study with the CTA 
will be the very few more distant objects with X-ray jets that are 
extended on several-arcmin - and thus typically 100-kpc - scales. 
The best example that we are aware of is NGC 625 1 (2 = 0.0243), 
whose large-scale X-ray jet, starting around 3 arcmin from the nu- 
cleus, was first detected with ROSAT (Mack, Kerp & Klein 1997) 
and later studied in detail with Chandra and XMM (Evans et al. 
2005). As discussed in the latter paper, it is not even clear whether 
the X-ray emission from the 100-kpc scale jet in NGC 625 1 is syn- 
chrotron in origin, and we do not have a detailed model for the 
electron spectrum or direct constraints on the jet speed or angle to 
the line of sight. However, we can make some rough estimates of 
the detectability of VHE 7-rays. Using the one-zone synchrotron 
model presented by Evans et al, we modelled region 1 of the ex- 
tended jet as a cylinder with a linear radial dependence of electron 
density and with 6 = 45° and /3 = 0.7, placed it at the appropriate 
distance from the host galaxy (which we assume for simplicity to 
have the same SED as M87) and computed the inverse-Compton 
emission in the VHE 7-ray band. We find a 1-TeV flux density at 
equipartition of ~ 7 x 10 -18 Jy, a factor ~ 20 below that measured 
by HESS for Cen A, and compute a predicted photon index of 2.5. 

3 There is some evidence that on large scales the energy density in FRI jets 
is dominated by non-radiating particles — e.g. Croston et al. (2008) — and 
so it is not completely implausible that the true magnetic field is greater 
than the field derived by assuming equipartition with the radiating electrons 
only, but without these CTA observations there will be no way of testing 
this model by direct observation. 



The CMB is by a large factor the dominant photon field here - light 
from the host galaxy, SSC and the EBL all come in an order of mag- 
nitude or more lower - and what prevents the TeV emission being 
much fainter as a consequence is the rather small break inferred by 
Evans et al. (2005) in the electron energy spectrum for the jet, cou- 
pled with its large volume, which means that the bulk of the jet's 
luminosity on these scales actually emerges in inverse-Compton in 
the 7-ray regime. Although the flux density we estimate here is 
only one value from the large range that we would obtain if we had 
the data to carry out more detailed modelling, the fact that it is only 
about an order of magnitude below what is already measurable for 
Cen A is encouraging, and suggest that NGC 625 1 and any similar 
objects might provide interesting targets for long exposures with 
the CTA. 



7 CONCLUSIONS 

The main conclusions from this paper can be summarized as fol- 
lows: 

(i) We have developed a framework for carrying out detailed 
inverse-Compton modelling of the kpc-scale jets of radio galax- 
ies in the regime where the anisotropic nature of inverse-Compton 
emission and Klein-Nishina effects are important, and applied it to 
the well-studied jets of Cen A and M87. 

(ii) We have shown that the predicted inverse-Compton flux 
density in VHE 7-rays on plausible models of Cen A is quite com- 
parable to what is observed. This means that the existing TeV ob- 
servations of Cen A put quite strong constraints on the magnetic 
field strength in the kpc-scale jet, particularly given that we have 
not modelled either the TeV electrons present in the compact knots 
in Cen A or the photon field from the presumed hidden blazar in 
its nucleus: B must be comparable to or larger than the equiparti- 
tion field, unless our beliefs about the probable jet speed or angle 
to the line of sight are very wrong. This is the first time that a TeV 
detection of an individual source has been used to put constraints 
on the jet magnetic field strength; our work supports the analysis of 
Stawarz et al. (2006b), which was based on the 7-ray background, 
in suggesting that the magnetic field in these jets cannot be much 
less than the equipartition value. 

(iii) Similar but weaker constraints can be derived for M87: it 
seems likely from our analysis that the large-scale jet contributes 
only a small fraction of the observed TeV emission from this 
source, which is consistent with the fact that the TeV emission is 
seen to be strongly variable. More observational work is needed in 
this case to pin down a non-variable component which could be 
compared with our predictions. 
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(iv) The prospects for study of extended X-ray jets with future, 
more sensitive VHE telescopes such as the CTA are limited primar- 
ily by the spatial resolution that will be achievable - in many cases 
kpc-scale jets may be detectable but unresolved from emission re- 
lated to the AGN or to the pc-scale jet. In Cen A, there is a very 
realistic possibility of resolving the jet, and we have shown that the 
possible 100-kpc-scale X-ray synchrotron jet in NGC6251, which 
should also be resolved by the CTA, may not be completely out 
of reach for such instruments. Angular resolution, rather than raw 
sensitivity, may be the limiting factor for the CTA in the field of 
extragalactic jets. 
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